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Abstract 


An  expans ic n  is  Given  far  the  5rcu11d.-r.tate  enemy  of  c.  cyst e:.  of 
interacting  I: cm  particles,  using  Ray  la  ig*- •  t -hrodl '  gc r  portuouuic n  theory 
with  the  scattering  part  of  the  interaction  ac  the  pertubrtion.  The  result 
in  second  crc.er  reproduces  the  energy  shift  found  by  various  authors.  The 
energy  spectrum  of  the  zero-order  Hard,  lucn  inn  c . n tains  a  gap  rf  nagnitude 
^0VQ  where  is  the  Fourier  transform  of  the  two -particle  interaction. 

This  gap  renoves  the  third-order  divergence  resulting  when  the  entire 
interaction  is  used  as  the  pertubaticn.  The  structure  factor,  computed  to 
first  order,  is  in  qualitative  agreenent  with  cxpe^ih.ent  for  large  uonentu:. 
values.  For.  non-zero  ter'perature,  expressions  are  obtained  for  the  depict ior 
of  the  zero-nonentu  1  node  and  the  specific  heat,  using  the  zero-order 
Hamiltonian.  The  specific  heat  exhibits  a  disc  ;ntir.uity  at  the  critical 
point.  Quantitative  estivates  for  the  magnitude  of  the  specific  heat  at 
the  critical  point  yield  results  which  arc  in  fair  agree. tent  with  empirical 
values  for  liquid  helium 
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Introduction 


As  is  well  known,  liquid  helium  II  can  be  idealize!  as  a  systen  of 
interacting  B  se  particles.  The  nor- inter c sting  system  'ras  studied  by 
London^"  who  showed  that  below  a  "critical "  temperature,  c.  condensation 
phenomenon  occurs.  The  average  number  of  particles  in  the  zero  momentum 
node  is  of  the  sane  order  of  magnitude  as  the  total  particle  number.  Iii 
contrast,  the  average  number  of  particles  in  any  other  state  is  of  the  order 
of  unity.  The  tern  "Brse-Einstein  condensation"  is  used  to  denote  the 
special  role  cf  the  zero-nonentun  node. 

For  the  interacting  systen,  it  will  not  suffice  to  use  simple 

2  3 

perturbation  theory.  Brueckner  and  Sawada  and  also  Huang  and  Yang  have 
shown  that  the  ground  state  energy  per  particle  is  not  independent  of  the 
total  particle  number  N  in  third  order  and  beyond:  instead,  it  becomes 
infinite  in  the  limit  of  large  N. 

The  first  microscopic  theory  of  the  interacting  systam  is  due  to 

k 

Bogoliubov  .  Assuming  sufficiently  weak  interatomic  forces,  he  showed  that 
the  Hamiltonian  could  be  diagonalized  by  a  canonical  transformation.  The 
transformed  Hamiltonian  represents  a  set  of  non-interacting  quasi -particles 
which  obey  Bose  statistics.  For  large  momenta,  the  energy  spectrum  of  the 
quasi-particles  is  quadratic  and  is  identical  to  chat  cf  the  ideal  systen. 

But  for  low  momenta,  the  spectrum  is  linear  (phonons)  a  radical  departure 
from  the  free  particle  results.  The  spectrum  found  is  similar  tc  that 
proposed  earlier  by  Landau^ ,  based  on  thermodynamic  data. 

Two  approximations  are  necessary  in  the  Bogolinbov  theory.  First,  it  is 
assumed  that  condensate  in  still  present  when  the  interactions  are  taken 
into  account.  Thus,  for  sufficient, ly  low  tempera tiucs,  the  average  value  of 
N  =  an+  a^  will  bo  large,  compared  to  unity.  (Ue  denote  the  number  operate" 
for  the  zero  mouenl’in  mode  by  N  ;  a  *  and  a  are  the  creation  and  destructi  ■> 

*  <  •>  y  »  >  ' 


operators  for  the  zero  n  orient  urn  node,  respectively. )  The  commutator 
£o>  »o+]  ,  (which  is  exactly  unity  for  Bose  statistics)  is  snail  compared 
to  .  Hence,  p.q  approximately  commutes  with  aQ+.  Both  these  operators 

can  therefore  be  replaced  by  the  constant  ^  N  'S^^;  the  value  of  ^ 
can  be  determined  later  by  the  requirements  of  self-consistency,  when  the 
ground  state  has  been  determined.  This  approximation  is  called  the  c -number 
approximation.  It  is  worthwhile  to  remark  that  the  c-number  approximation 
ignores  the  fluctuation  in  N  =  a  +  a  ;  this  is  justified  because  the  theory 

O  CO 

yields  a  fluctuation  in  Nq  which  is  of  the  order  of  magnitude  of  '  , 

much  smaller  than  ^  • 

The  second  approximation  is  to  omit  the  terras  in  the  Hamiltonian  H  which 

+ 


contain  less  than  two  factors  of  a  cr  a 

u  o 


A  straightforward  analysis 

shows  that  this  approximation  is  valid  only  for  weak  interactions.  That  this 

is  to  be  expected  can  be  seen  from  the  following  argunent:  Due  to  the 

interaction,  states  other  than  the  zero -momentum  node  are  occupied.  It  is 

said  that  the  zero-ncnentun  node  is  "depleted".  For  sufficiently  weak 

interactions,  the  depletion  is  small.  Then  the  terms  of  H  containing  less 

than  two  factors  of  aQ  or  a^+  will  be  unimportant  since 

Lee,  Huang  and  Yang^  have  investigated  the  hard  sphere  Bose  gas,  using  a 

pseudopotential  for  the  two-particle  interaction.  Their  method  can  be  shown 

to  be  equivalent  to  that  of  Bogolinbov;  their  results  for  the  ground  state 

energy  and  the  excitation  spectrum  are  similar  to  the  Bogolinbov  results.  Lee 

7 

and  Yang  have  extended  the  study  to  nen-zero  temperature  and  have  computed 
the  specific  heat.  The  specific  heat  is  found  to  exhibit  a  discontinuity  at 
the  critical  temperature;  the  magnitude  of  the  discontinuity  is 
proportional  to  the  hard  sphere  diameter. 


Extensive  studies  have  been  completed  on  the  interaction  between  the 
quasi-particles.  (The  quasi -particle  interaction  is  contained  in  the  terns 
dropped  by  the  Bogoliubov  truncation  of  the  Hamiltonian. )  The  ground  state 
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energy  has  been  examined  by  Wu  for  the  herd  sphere  gas.  Applying 
pertubation  theory  to  the  quasi-particle  interaction,  Wu  obtained  the  next 
ordor  tern  beyond  that  calculated  by  Lee,  Huang  and  Yang. 

9 

By  euplcying  the  Green  function  method  of  quantum  field  theory,  Beliacv 

and  later  Hugenholtz  and  Pines'^1  have  confirmed  the  results  of  Wu. 

Hohenberg  has  reformulated  the  Green  function  method  using  variational 

derivative  techniques.  His  work  enables  the  computation  of  the  denisty 

correlation  function  to  one  order  bey  nd  the  Bogcliubov  approximation. 

Several  workers  have  studied  the  properties  of  truncated  Hamiltonians 

12 

which  are  more  complete  than  that  of  Bcgclinbov.  Girardeau  and  Arnawitt 
have  added  all  pair  scattering  terms  to  ‘the  Bcgoliubov  Hamiltonian  and 
diagonalized  the  resulting  Hamiltonian  by  a  canonical  transformation.  The 
energy  spectrum  that  is  obtained  contains  a  gap  in  the  limit  of  long  wave¬ 
length  and  thus  shows  that  the  sucess  of  Bcgoliubov  in  obtaining  a  phonon 
spectrum  was  fortuitous. 

13 

Takano  J  has  considered  correlations  beyond  simple  pair  scattering.  He 
shews  the  existence  of  a  collective  mode  which  removes  the  gap.  H  wever,  it 
is  not  easy  to  obtain  the  corrections  to  the  equation-cf-noticn  used  by 
Takano. 

14 

Luban  has  considered  the  thermodynamics  of  the  Girardeau-Arnowitt  model. 
The  phase  transition  at  the  critical  temperature  is  found  to  be  of  first 
order. 

Lieb  end  Leninger^  have  shown  that  in  one  dimension  a  system  of  bosons 

interacting  pairwise  via  delta-function  potentials  can  be  solved  exactly. 

16 

Lieb  has  investigated  the  excitation  spectrum  and  has  shown  that  in  this 
model  there  is  a  double  spectrum.  The  two  dispersion  curves  possess  a 
common  slope  in  the  limit  cf  low  momentum. 

The  work  described  in  this  report  has  been  motivated  by  the  desire  to 

obtain  a  systematic  perturbation  theory  to  compute  basic  properties  of  the 
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Bose  systen  in  a  region  of  temperature  above  as  well  as  below  the  critical 
temperature.  For  this  purpose,  we  develop  a  method  which  avoids  the  use  of 
the  Bogolinbov  d- number  approximation. 

Section  I  begins  with  a  discussion  of  Rayleigh-Schrodinger  perturbation 
theory  (R-S)  applied  to  the  Bose  gas.  We  show  how  the  straight-forward  use 
of  R-S  perturbation  theory  in  a  calculation  of  the  ground  state  energy  leads 
to  divergent  results  in  third  and  higher  orders.  A  simple  modification  of 
the  R-S  method  is  given  which  avoids  these  divergent  tens.  The  structure 
factor  is  computed  to  second  order  in  the  modified  perturbation  scheme. 

The  calculation  reveals  the  expansion  parameter,  which  is  a  function  of  the 
Fourier  transform  of  the  potential.  The  parameter  is  small  at  large  momentum 
values  but  approaches  one  half  in  the  limit  of  vanishing  momentum.  The 
structure  factor  is  in  qualitative  agreement  with  experiment  in  the  region 
of  large  momentum.  In  the  low  momentum  region,  the  perturbation  series  for 
the  structure  factor  converges  slowly  so  that  the  first  few  terms  do  not 
provide  a  good  representation. 

In  Section  II,  we  present  a  mathematical  introduction  to  the  methods  used 
in  our  calculation  of  the  statistical  mechanics  of  the  lowest  order 
Hamiltonian.  The  ideal  gas  is  used  as  an  illustration. 

Section  IH  contains  the  calculation  of  the  statistical  mechanics  of  the 

Bose  gas  interacting  through  the  Hartree-Fock  potential.  A  discontinuity  in 

the  specific  heat  is  obtained.  The  magnitude  of  the  discontinuity  is 

proportional  to  the  rate  of  depletion  of  the  zero-momentum  node  with  increasing 

17 

tenperature.  By  using  the  Yntena-Schneider  potential,  qualitative  agreement 
with  the  experimental  specific  heat  is  obtained.  The  range  of  validity  of 
the  calculation  is  examined  and  is  found  to  exclude  only  the  case  of  a  weak 
repulsive  interaction  near  the  transition  tenperature. 
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Section  I 


Modified  Fcrtubation  Theory  and  Structure 
Factor  of  a  Many  -  Boson  System 


Ue  first  wish  to  study  the  properties  of  a  system  of  pj  int-r  ''tin,'-; 
Bose  particles  with  the  use  of  Rayleigh  -  Sehred'n ^cr  (R-S)  periuvt.'": ' '' " 
theory. 

The  Hamiltonian  we  shall  consider  is'' 

H  =  H  +  H 


1C 


i  ? 

H  -  I  7-  0.^0 
°  A  *  A  4 


H  -  -U  J  V  nr 

'  ^  r;  w  V n  Q-n' 


t 


Cl, 


(1.1) 


a  o 


M,*'  *  ^  * 


In 


Eq.  (1.1),  T  =  T,3- A1/  3,^1 


is  the  kinetic  energy  cf 


particles  of  momentum  j)\:  and  mass  m.  The  operators  a*  ana  are 
creation  and  annihilation  operators  for  these  particles  and  satisfy  the 
commutation  relations: 


C  a  a  f  1  -  f 
r  5  *  J 


C  a  ,  or  ]  =  *  • 
A‘  >  A  > 

<1.2) 


V  is  the  Fourier  transform  of  the  two-body  potential  V  (  Y  )  J 


1r 


\/  -  <? 
*  J 


The  allowed  values  of  the  wave  vector  in  the  suras  are  determined  by 


imposing  periodic  boundary  conditions  in  a  cube  of  length  L  end  volume 
■R-  =  L3.  Hence,  Jk  ~  2  IT  W /{_  A  :=.  ~1  TT  /V\ ,  1 1 


’V  i7r/H,i/<- 


/TA  -r»7 


syy\ 


~  O  st  I 


•s 


In  calculating  the  ground  state  energy  of  the  Hamiltonian  II,  a 
integral  appears  in  the  third  order  of  the  perturbation  expans. A'.  >* 

section  l.A  we  show  how  this  difficulty  can  be  circumvented  by  c  iificaticn 

o*’  the  R-S  expansion,  and  we  give  results  for  the  ground  state  e  y  up  to 

third  order.  In  section  l.B,  the  structure  factor  of  the  Bose  s v~to.  is 
computed  in  the  modified  perturbation  scheme. 

A.  Modification  of  the  R-S  Perturbation  Expansion 


The  interaction  (Eq.  l.l)  contains  two  types  of  terras.  First,  we  hnv. 

->  > 

terns  that  are  called  forward  ar.d  exchange  scattering  terms.  They  correspond 
to  £  =  0  and  £  =  k’-k.  These  terras  do  not  result  in  any  net  exchange  of 
momentum  between  the  Bose  particles.  Hence,  the  plane  wave  eigenstates 
of  Hq  are  also  eigenstates  of  the  forward  and  exchange  terras.  We  are  thus 
led  to  separate  out  these  terras  from  the  interaction  H^  and  include  then  ir. 
the  unperturbed  Hamiltonian. 

The  second  type  of  term  in  H^  is  composed  of  that  part  of  the  interaction 
which  remains,  after  the  forward  and  exchange  terms  are  removed.  These  terms 
allow  the  Bose  particles  to  exchange  momentum.  These  terras  will  be  the 
perturbation  and  are  denoted  by  H^' . 
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H,  = 


Sl 


where 


i  z  , 

/  A  A 

l*° 


V  + 
-i  V 


a.  <Z  q 

'?  *  A'  I 


V  - 


a£  q-A  1 


Examinf ng  the  modified  free  Hamiltonian  H  '  v.?  find  that  the  ;■/ 

o 

energy  is  not  zero,  but  is  given  by 


1  u 


E  '  -  H  OH)  -t/ 

o  V 

p_  ° 


Eq.  (l .4)  is  easily  obtained  by  simply  nettiv 


v;  n 


N  ar.d  - 


\  ^  O  The  next  excited  state  of  the  Hamiltonian  H  '  with  total 

c 

/  — N,  — )  V 

momentum  zero  is  the  state  cf  a  single  excited  pair  with  momentum  (q,  -  qj. 
The  energy  of  this  state  is  obtained  by  setting  r%  =  Af  ^ 


m  -  j 
%  > 


rr\  =1  ) 

~c.<  -> 


rnfi.~  °  %,-T- 


Thus, 


(K). 


f 


it 


4- 


lj  OV-Qy  +  +  \fz, 

<t  J),  ->  4j 


sz 


(1.5) 


If  we  take  the  difference  between  (.1.5)  and  (1.4),  we  have 

O  =  iil  +  iiELijy  .  (.,.5) 

In  the  limit  q  — )  0,  then'  Jr-  a  gap  of  m. '.gait  ado  ;l(/V*2)  /  t.  betv 

the  ground  state  and  the  first  excited  stale  of  the  Hamiltonian  H  1 .  It  is 


the  presence  of  thir.  gap  in  t.he  energy  for  an  excited  pair  that  nrevents  'he 


divergence  of  the  third  order  perturbation  term,  as  ve  shall  soon  see. 


Using  Hq'  as  the  unperturbed  Hamiltonian  and  H^'  as  the  interaction, 
the  energy  corrections  to  the  ground  state  are  as  follows: 


£#  =  fS(M-/)y 

°  a  jx  0 


£ 


-  o 


£  '  -  M  (N-i) 

a  SI  C^rrf 


; 


V  4 


%■ 


o 


N 


j_  +• 

'XT 


2  (/y-z')  \r  J 
J'L  '& 


g  /  4- s  got,  .eg 


=  rtW 

"  '  £  -f  IM-ilVi/*-)  C*'V~, 

Note  that  the  denominations  in  the  integrands  are  cf  the  form 

^Jr1  +- 

<3.  f  1  & 

possibility  of  a  lev  r.oncn+u.i  div.-.rg' nee. 


Tne  second  torn  precludes  the 


i  B.  Structure  Factor  cf  Interacting  Bose  Gas 

\ 

25  25 

It  has  been  shovn  by  Van  Hove^  and  also  by  Cohen  and  Feynman*-  that  in 
an  inelastic  neutron  scattering  experiment  one  can  r .ensure  the  probability 
per  unit  time  that  a  slow  neutron  transfer  energy  ~t)us  and  Momentum  to 

the  Boson  system.  For  the  ground  state,  the  probability  is 

-  ft  5  68,  ur)  (1'8) 

in  the  Born  approximation.  Here,  A  it;  a  constant  which  characterises  the 
neutron  -  Boson  interact i-.n,  whllo  !'<  (  &  m  )  is  cal  led  the  dynamic  structure 
factor  and  depends  only  nit  the  prnpm  i.  !ou  nf  the  l1  ::nu  :iyr.u_r.:.  f:  is  define1 


(1-9) 


=  2  (/)  U^-urJ  . 


/  -A  /no 


In  Eq.  (1.9), 


t 


?A 


T 


is  the  density  fluctuation  of  momentum  k  ,  de.r 


fA  '  $  d\ 


^  -q 

<  A-v  5  ^ 
=  2  C" 

j 

<r 


'-A  ~j 

<  -A  y 


o 


KT  "0 

-  >  Cl  Q 

f  * 

t 


(i.io) 


The  subscript  n  and  0  in  /  P  )  denote  matrix  eleuer.ts  between  the  exo.^ 

(  /a  Vic 

boson  state  n  and  the  ground  state.  Also,  ~t)  uC»c  is  tuo  cu.  Ltr.m'cn  ener.- 
of  the  state  n  above  the  ground  state. 

The  differential  cross  section  for  inelastic  neutron  scattering  can  be 
obtained  fron  Eq.  (1.3)  by  integrating  ever  the  energy  Thus,  the 

differential  cross  section  is  prop'-nti'anl  t  ' 


5  i,  /  V  S(A, 


u/ 


(1.11) 


o 


J0„  is  called  the  Diquid  structure  factor.  It  ur  the  Fourier  oransf om 


the  pair  correlation  functicn  .Jb  (  P  )•  Thus, 

0 


(k 


}  J  Jiff,  ffi 


it  Jl  ff(p)  fct)  lu\i 


-<a  v 
) 

(1 .1') 
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The  state 


refers  to  the  exact  ground  state. 

To  obtain  Sd ft  we  use  the  perturbation  expansion  of  the  ground  state 
for  the  natrix  element  of  the  density  fluctuations  in  Eq.  (1.12). 

Frcn  R-S  perturbation  theory,  the  ground  state  can  be  expanded  to  second 
order . 


The  result  is 

I 

l°> 


<M°>  + 


>  I  oy 

*  £  '  -  £  ' 


-+  Z  l^><  ri  '7Y\  /«>> 


/yv) 

) 


.U-iM 


•'*  efy>  / 

The  coefficient  a  of  the  unpertnrbatcd  state  /c)  in  Eq.  (l.  14)  is  c. 


normalization  fact  <r  which  .must  be  chosen  such  the. 
This  calculation  yields 


a  =  -i 


z  i/V;j 


* 


^7 


a 


Eq.  (l. lU)  can  be  applied  specifically  to  cur  Hamiltonian,  Eq.  (1.3). 

-f-  + 


A  straight  forward  calculati-  n  yields  the  result 

U1  -  / 


°y  —  a  |o^>  -  y 


aa  £ 


%  /  Tv?  -h  2.  jy 


^  lA0v-j)(/Y~2 )^1~3) 
4  J7.a 


_ ^rfr  t  -f  ■ 

w  K  a  <1  a 

Z_  L/_ 


t  ~h  -t  t 

a  10) 

Hr  '1 


i'.fr  ][£+£/} 

1  I  -*n  1  * 


J  J 


(The  notation  j  ^  refers  to  the  ground  state  of  the  non- 

interacting  system  of  N-2  particles,  etc. ) 

t 


To  evaluate 


S(A) 

t 


(V  H 


l°>, 


ve  use  the  identity  p 

J  A  'A 


^  +  2  af ,  a  a  a 

JJ"  J-A  4'  jp  ^ 


an  immediate  consequence  r  f  the  corxiutation  relations  Eq.  (1.2). 

Thus  .  /  \  c  A  j 

su  ?  ,  t  +  s  H<°l  VV°)h- 


n  1  Jr  *-£  J+a  J}‘ 

u? 


i  y  ;  _ 

Inserting  cur  result  f<  r  |  O  ^  in  the  at  ve  expression  for  -->(£_) 
yields  cur  result  for  the  structure  factor: 


suU)  "  i- i  / 


J  ^  k\  u 


■+ 1  />  X 


A 


a  J 


v 


yt-j  +$%,K  _ 


Examining  Eq.  (l.l5),  it  is  seen  that  the  expansion  parameter  for 
can  be  roughly  designated  by  —  P  w?  J—  ^  f  P  Hence,  at  high  k,  a 
rapid  convergence  ■ f  the  series  can  be  expected.  However,  at  email  k,  the 
series  converges  slowly,  since  each  tern  in  the  series  is  r<  uglily  one  half 
the  preceding  terv.  in  the  limit  k  —)  0. 

Suppose  we  retain  only  the  first  order  correction  to  S ($)  ir»  Ecl-  ( -1- •  1 5 ) 
and  compose  the  qualit.nl,  j  ve  t'om  of  the  resu1  i  s  with  experiment.  The 
experiment.''.  1  values  for  S(f)  are  ciiarac ter ixed  by  damped  nsc  i.1  Let  L,  : ab 
unity  for  high  values  f  h.  Experimental  curves  f.  r  S/A)  aro  contained  in 


27. 

the  W'  rk  performed  uy  flenshaw  . 
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To  compare  with  Eq.  (1.5),  suppose  we  employ  the  soft  core  potential 

~  ^  -A  <l 

(aT)^ 


-  (Lera. 


■A<q 


where  a  is  the  radius  cf  the  soft  core  and  where  ~\J  -  C  Q.  J ^  ^ 

(Here,  C  is  the  magnitude  of  the  potential  barrier. )  Then,  fer  large  It, 
Eq.  (l.l5)  becomes 


(*) 

S  (A)  -  / 


t  ~Aa~ 


ruk 


— 

-  A  Q 


(1.16) 


Note  that  the  general  form  of  Eq.  (l.lo)  is  in  accord  with  the  empirical 
results . 

In  conclusion,  ve  night  add  that  we  expect  the  series  expansion  for  the 
ground  state  energy,  Eq.  ( 1 . T )  to  be  rapidly  converging.  For  the  phase  space 
factors  in  the  integrals  over  the  momentum  q  in  the  expressions  for  the 
correct! >ns  to  the  ground  stace  energy  are  such  as  to  deemphasize  the  low 
momentum  region. 

So  the  slew  convergence  of  the  series  for  S(A)  >  in  the  region  cf 
small  k,  has  no  counterpart  in  the  series  for  the  ground  state  energy  because 
in  the  latter  series,  the  sun  over  all  momentum  levels  contributes  a  density 
cf  states  factor  that  diminishes  the  importance  of  the  lew  momentum  states. 
Hence,  we  believe  that  the  expressions  Eq.  ( 1 . T )  for  the  ground  state  energy 
up  to  second  order  is  a  good  approximation  to  the  exact  perturbed  ground 
state  energy 

In  the  next  section,  we  shall  pursue  an  investigation  cf  the  thermo¬ 
dynamic  properties  of  a  Bose  system  using  the  lowest  order  Hamiltonian  H', 
with  the  view  of  obtaining  reasonable  quantitative  agreements  with  experiment 
for  quantities  such  as  the  specific  heat. 

-Ill- 


Section  II 


Base  Gas  with  Hartree-Fock  Interaction 
A.  The  Transition  Temperature  and  Zero-Momentum  Mode  Occupation 


According  to  statistical  mechanics,  all  thermodynamic 
properties  can  be  obtained  from  the  canonical  partition 
function  (O)  ,  defined  by 


'I 


Here,  cr.»,  is  the  th  energy  level  of  the  K-body 
system,  T  is  the  absolute  tempera  cure  and  K  is  Boltzman's 
constant. 


A  partial  list  of  thermodynamic  quantities  of  interest, 
derivable  ficr  Q  ,  are  giver.  be  lev/: 


ft 


(HeL"iholtz  Free 
Energy ) 


V 


>3 


Z-<>. 


1 


(internal  Energy) 


c 


v 


) 

<?r  4  l 


(Heat  Capacity  at 
Constant  Volume ) 


To  bee in  the  tack  of  ev 
we  insert  the  eigenvalues  of 
for  the  values  of  £  in 

->i 


H'a 


the  partition  fuiiction 
•  Gtvon  in  F fr  H  3) 
Thus.  0 


*\zfl 


i^a- 


Each  state  in  the  canonical  ensemble  has  a  fixed  total 
particle  nui  :ber  /\|  .  Hence  the  second  term  in  the  exponent 

can  be  removed  outside  the  sujfl.  The  partition  function  can  no;/ 
be  written 


-  (J  tyn)- 1?  VQ 


-JU  n-  /  ) 


(2.3  ) 


where 


(y)  is  defined  by 


$2  - 


'T  -f-  2 


r  a 

^  1  ?  ' ) 


To  cor 


s' 


t->’<  S  * 

ve  'r.  .".no  (j/  •'  .  |  'ft  )  by 


IL'V/j?-  n; 


A! 


A 

A  i  o 

5  - 


-i  s  v  V'O} 

*1* 

N 


wnere 


A'! 


Hence, 


QM  -  r  9 


£  ^ 


^1=0 
J 


—  1 5 1~ 


Note  that  can  be  thought  of  simply  as  the  partition 

function  of  a  system  of  N  particle,  whose  single  particle  levels  do  not  include 


the  zero  momentum  level 


a) 


M  A  X 


denote  the  m  ri'  urn  term  in  the  cum 


on  ft  in  Eq.  (2.6).  Then,  in  the  limit  of  large  N 

O 


(2.7) 


The  verification  of  Eq.  v1.7;  is  simple: 


Write 


Thus , 


^1;  ^rwflX  7  -"rCj— 7-' 
N - -  « 


"^7  77  rtf  +  /f%H 

~A  N  N 


M 

o 


is  not  ser s  in  the  limit  N  5  02  f  hue 


J  Q/j  —  Ly M fix 


which  verifies  Eq.  (2.7). 


S  /■  \ 

The  generating  function  for  Q  (  ft  zft  y  ic  defined  by 


=2  fe,  ^  =•  2 


rf\ ,  ^ 


(2.8) 


(  ?  0  =  7- 


/  y  ^ 


•y  >1  2* 

i  M  i 


#  V  J ■^■■o  ™$j  'ef  J 

'i  -Kk 


The  prime  on  the  summnt  i<  .n  sign:;  i ml  Inil  < •:;  i  h.  ili-lel. ton  of  the  zero  momentum 


.-.-ado  from  the  Tin'  yml.t.i  ^  'Tt  J  d.-ii.>te..  tin-  f.  -t.  that  we  sum  over  the 

entire  set  of  oecupat  u  ,.i  a  mil.. -rs  (<  | « i  /;  ),  cuch  mvupnt  ion  number  assumes 


the  values:  1,2,  i 


,  O 


To  evaluate  Eq.  (2.9),  we  replace  '??  by  <^0'!  /  ,  vhere  the  average 

%  v 

is  performed  in  the  ideal  ;;as  ensemble .  and  will  be  computed  later.  This 
replacement  enables  the  calculation  of  ^he  s  ims  in  Eg.  (2-9)  since  the 
summands  are  now  separable.  We  obtain 


(2.10) 


IIov  can  be  obtained  from  (.'2:')  by  means  of  the  identity 

~  2F4  o  —^7  ^  (2.11) 


an  immediate  consequence  cf  Eq.  (2.3).  The  contour  is  any  closed  curve 
containing  the  pole  •/,  ~  0  and  lying  within  the  radius  of  convergence  of 

2  ^>.  ^ 


C?\ 


In  the  limit  N  — ^  ^  y 

Eq.  (2.1)  can  be  written  in  the  fom 

R  F  1/h  *) 


Qm  "  IT,-  §><&  & 


(2 .22) 


where  F  C^^')  ~  ~f  ( f)  ^  ^  >P  '^3  ^ 

-f  (  P  ?  )  —  ^ 

^  ’  /V  c«,  ^ 

The  integral  (h.Jf)  I:;  in  a  r.  >rm  in  which  il  can  1h  evaluated  by  the 
Method  of  saddle -(joint,  i  n(  •  •gr.ot.  Imi.  (iv-talla  can  he  feiiini  in  reference  23, 
particularly  pager.  >1 M  In  5.  ni:...  1  .-:<  i  <  n.  e  .'>»  Chapter  1.  for  an  outline 


l  i 


of  the  general  ineUio.l). 


The  result  is  that  is  giver,  "by 


JL 


N  -  s-i<y-q  P{ 

zh7L 


/  /V  (_  Or)  _  Jy  JL'-'j  7* 


./ 


(-.13) 


where  z  is  the  solution  of  the  equation 


2 

<?  ir 


Cl ) 


-J~L 


=  f 


(2.14) 


Substituting  Eq.  (2.10)  in  Eqs.  (2.13)  and  (2.14)  gives  the  result 


g  ch-*qi%')  _  (*, 0  (2 


iM 


where  z  is  determined  from  the  solution  of  the  equation 


f>o-s)  =  l t 

o?rr' 


(2.15) 


We  now  wish  to  naxl.dze 


Q(H 


3-/\ 


■oo,,  /yi 


J 


to  obtain  from 


Eq.  (2.7).  The  maximal  condition,  using  Eqs.  (2.l4)  and  (2.15)  is 


tf5  *  </V+h  5  w titu: 


Tire  two  simultaneous  eqirations  (2.  1  3)  and  (2.15)  determine  z  and  /  . 


Eq.  ( 2 .  l^i )  then  gives  0^. 


We  must  now  evaluate  the  integral!-..  First,  consider  '.he  sura 

T’o  a  ^  v  y  which  appears  in  -ho  exponent  of  the 

L  rL  1 A  #  o 

integral::  in  Eq.  (2.15)  and  (t’.'h).  Tin-  region  in  which  q  large  contribute. 


negligibly,  due  to  the  exponential  factor.  We  therefore  replace  ^  by  ^  , 


Changing  the  sun  to  an  integral  yields  the  result 

i  z  V 

fa  7-*  J 


i  Z  V,  <*.>  --  U  «■  /r, 


<r  = 


V/T 


) 


Vi 


oo 


F  (*)  =  n 

"Vi  / 


cJ~X 


0 


(2.1?) 

'T?  “  / 


r 


«> 


-)  '  -H  I  W 

The  numerical  value  of  F,  ?  I  )  =  2.,  6 /Z  . 

'L 

Defining  by  ^  q~  ^  we  can  now  rewrite  Eq.  (2.15) 

and  Eq.  (2.16)  as  follows: 


; 

oo  i 


po-s)  =  j £j£l _ 

<D  ^)t  $  jfVQ  + '  S?*  |  ^  (2.13) 


-L 


-  -- 


^‘h 


(2.19) 


The  chemical  potential  Z/  is  defined  by  ~  -  °^/(3  • 

To  proceed  further,  we  define  the  temperature  T5  by  T°  =  Lf  TT  h 


,  V3 


^  A 


f 


VJe  also  define  the  dimensionless  parameter  x  by 


and  the  parameter 


The  magnitude  of  T  for  helium  is  5.9?°K. 

r  oo 


-  A/ 


/T  I 


y-  -.21 


x  (1 /To)'U 


Wo  also  define  F  (s)  = 

n  '  ' 


Pfr  ) 


3  ^ /(e  r*-i) 


With  this  notation,  Eqr .  (2.18)  and  (2.19)  become 


^YV i-o-  x 


T 


) 
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(2.20) 


and 


*  =  ( ^r-rx(L<l\r,-*)) 

f  ’  K  ' 


Of 


(2.21) 


can  be  eliminated  from  the  pair  of  equations,  Eq.  \2.20)  and  (2.2l). 


Tlie  result  is 

ry* 

r  ' 


(!-<>)  -  F  fj/xT 


!k 


r+%- 


}~  x  7 


(2.22) 


Eq.  (2.22)  is  a  transcendental  equation  for  J  the  fractional  occupation  of  cue 

zero  momentum  made.  For  each  value  of  x,  there  will  be  a  solution  for  $ .  Of 

course,  x  —  P  Vo  is  a  measure  of  the  strength  of  the  interaction. 

»r7" 

Tr.c  transition  temperature  T^  is  defined  to  be  the  value  of  T  for  which 
^  =  0.  Hence,  from  Eq.  (2.22),  we  have 


T_ 

"T 


=  F 


-J-  \U  -  T “O  \ 


-  A  / 

r. 


as  the  equation  giving  . 


We  solve  Eq.  (3-25)  graphically  for  ~^~<z  jy  as  a  function  of  x,  by 

imparting  values  to  t  vr  and  determining  the  corresponding  value  of 

e 

x.  (Numerical  values  of  wcre  obtained  from  reference  1,  page  U4.  ,  The 

result  is  a  single-valued  function  of  x.  which  we  iJlustrate  in  Fig.  3-1. 


/  o 


(2.23) 


; 
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- — - 

As  x  approaches  infinity  in  the  positive  direction,  j ~J~  ° 

approaches  the  value  0.8^  fro::!  above .  As  x  approaches  infinity  in  the  negative 
direction  /c  j  ~f  up,  reaches  the  value  0.84  fron  be  lor.  Ac  x  .approaches  zero 

from  negative  values,  ~T~r  J~~T'  CL^  $2  7  “c  ^hat  in  this  case,  approaches 

the  ideal  gas  transition  temperature  ,  .$27  Ta.  ITote  that  the  effect  of  the 
interaction  is  always  to  raise  the  critical  temperature. 

For  O  <_  X  .  8  6  '~7~  /  7~  is  greater  than  two  and 

?  C  /  '  o 

rapidly  increases  with  an  almost  vertical  slope.  We  shall  show  presently  that 
the  values  Tc  /T°  >5,  (corresponding  to  O  (  .  i 6  )  are  net 

bonafide  critical  tempo rat urei and  will  also  discuss  the  criteria  for  an 
ac.  eptable  critical  temperature. 

The  density  of  liquid  helium  in  the  \-regior.  is  of  the  order  of 


„21 


.10  ‘  cm.  ,  and 


I  for  this  density  AT T  —  /o 


~t  6 


ergs.  Inserting  these 


vaxues  o. 


P 


and 


AT  T°  into  the  definition  of  gives  x 


10 


3 


(In 


the  calculation.  we  have  used  y/^  ~  / O  fa  ,  a  result  obtained  from 


numerical  estimates  based  on  t>ho  Yntema-Schneider  potential’*’1  ). 

We  now  examine  temperatures  below  /  ~  .  In  Eq.  (2.22),  wo  can  expand 

v_ 

1 


the  function  F 


3/r 


for  small  values  of  the  argument.  The  expansion  is 


?-n  ~  o  srv\  I 


where  b  is  the  Riemann  zotu-funcLion. 

o 


Employ  in,"  this  expansion  (which  converges  rapidly  foi  small  £  ) ,  we  obtain 


In  obtaining  Eq. 
Per  x  ^  0.85, 

Then  Eq.  (2.22)  becomes 


(2.2?i ),  we  have  also  expanded  in  powers  of 
the  large  argument  expansion  of  ^  (°f ) 


{. 

1 

can  be  used  . 


(2.25) 


Before  turning  to  the  calculation  of  thermodynamic  quantities,  we  examine 

the  range  of  validity  of  our  results  for  fT  .  We  must  require  that  (  I  AA  ) 

C  o'T  c 

be  negative  for  a  true  transition.  (That  is,  5  must  be  a  decreasing  function 

of  T  as  T  — T  .  )  If  we  compute  C.T  <5/  )  we  find  that  this  requirement 
c  '  3  T  X 

does  hold  except  for  x  .85.  Thus,  the  results  are  valid  for  all  x  except 
near  the  transition  temperature.  Near  the  transition  temperature  the  results  are 
valid  for  negative  x  and  for  2  .85  but  not  for  x  .85. 


B.  Thermodynamic  Function.; 

Wc  are  now  in  a  soniHo:  i.<>  .v.  Ion  Ini  o  tin*  enor,;/  and  '.pecifie  heat  of  the 
system  of  Bose  particles  inienvt  lug  with  the  Hnrt.ree-Fock  part  of  the  interact!  i 
'Hie  internal  eiioury  >>l'  i.Ih-  sy  leu  I  ■  define.!  by 


U  = 


<2 


A  ^  1  •» 


(2.26) 


Using  Eqs .  (2.3),  (2.10)  and  (2.lU),  we  have 


%  = 
~^r 


-t/V 


r  AO 

1 M  i  ^  /  -  d'.y  ^  (^  -  //  £ 

+-i  T'7"  )+6-$)otJ> ' 

ZJL  £  -+j>.}  i~A  *  ' 


(2.27) 


Combining  Eqs.  (2.26),  (2.27)  and  using  Eqs.  (2.15);  (2. 1.8)  and  (2.19) 


:ives  the  result 


i  z  /V*  +1(ZY(ht)F 


7-t±  1zV<\>\+^) 
0  2n  /  '  ^  7 


•4-  •?  $  p  V  -+ 

/  O 


% 


F  X  <■"* > 


/O'/; 


(2.26) 


To  obtain  the  specific  heat,  ve  differentiate  Eq.  (2.28)  with  respect 
to  the  temperature,  mate  use  of  Eq.  (3-ll),  the  fact  that  J?  F ^  (c><  )  — 

t 


md  Eqs.  (2.20)  and  (2.17).  The  result  i. 


/y  n  '  ^  7  -V 2_  '  7“  ^  To/  7-  ^ 

+  20-).)^-r  5?  )  xT%  [W>\ >/T  f,  ,US) 

f  ^  T  I7O  Jhi 

+  K'-  ))  t  *7.°  ■/ jf/- ^  /  A/r  V'i  .  ,  _  x  7'  ^Ar  ~J  . 

-7 -  7  V  c‘  ' 

Eincc  7”  <7j  i:;  I  i  .:>•• 1 . «u.;  :.t  1  hr  \  ■print  .  \.V  IHV  l’lVI:.  Eq.  (2.29) 

c)T 

t.h'.l  the  specific:  heat,  I : .  •)  i:  ••.ml  li.n.  >u  ;  »•!  1 1 1 1  I  ■■  i|h  m .  mv. 


To  examine  the  result,  Eq.  (2.29) 
temperature .  An  r->  O  .  5  ~>  1 

]  '  —  r 

calculation  yields  the  resu1 L 


first  consider  the  limit  of  low 

A  len.jt.hiy,  hut  straight- f'  award 


C. 

V' 

1  Tn  /. 


r-^c 


(2.30) 


The  exponential  result  is  expected  since  there  is  an  energy  gap  in  the 
energy  spectrum. 

Next,  consider  the  care  in  which  the  temperature  is  in  the  neighborhood  of 
the  critical  temperature.  For  temperatures  above  T  .  I  -  &  i  =  C  and  1^1  - 

JT 


(2.31) 


At  the  critical  temperature,  is  essentially  a  constant,  independent 

of  x.  ( C y ^  j  /\{/f  at  Tc  varies  fren  a  value  of  3*5  at  x  1  to  a  value 

of  3-2  at  x  =  ex’ 


For  temperatures  just  below  T(>,  Eq.  (2. £9)  yields 


1 

.  A 


(2.1  ) 


We  can  now  compute  the  magnitude  of  the  diseonu'  unity  cf  the  specific  heat 
o.t  the  transition: 


C 

■V 

ITT 


r:  -T-o  _ 


j  (  «T/ 

(2.33) 


The  discontinuity,  Eq.  (2.33)  is  proportional  both  to  x  r.nd  to 


tTc  sketch 


cu  /hh 


in  Fig.  2  for  large  x. 


Conclus ion 

The  exterir.ental  valuer,  '’or  c  at  the  X-ncirt  have  been  determined  by 
1  v 


31  /  o  , 

Fair-bank,  Buckingham  and  Kelierr  (to  within  10  degree  of  -he  > -point ) .  Their 


results  indirate  that  c  is  finite  and  continuous  nt  T,  and  that  c  has  a 

v  X  p 

1  ugr.rithnic  singularity  at  T  . 

A. 

The  theoretical  result  for  / /\)  !\  J T  ’  liven  by' 

s  in  order  of  magnitude  agreement  with  tin  observed  value  n" 


] '•! •  (2.32)  i 


(for  pressure  of  0.05  atmosphere),  as  given  in  reference  31-  However, 

Eq.  (2.33)  predicts  n  discontinuity  in  (r  / EK)(  ,  disagreeing  with  the 

results  contained  in  ••r-forc-  ice  31-  But  the  discontinu  .us  behavior  of 

c  /JJK  at  the  cr:it  i».r  1  roirt  is  in  oualit^t  h/o  --  eero’ t  with  experi- 
v 

rental  data  on  c  / I-.1C  reportov..  in  /A  '-.iff.  .  hue  cru-e:  "mental  values  of 

v' 

32 

c  noted  by  Atkins  are  derived  f rc :  work  performed  prior  to  the  work 
v 

of  Fail-bank,  Bucking  bar  and  Keller-?  The  theoretical  expression  for 

(c  )  /lK  in  Eq..  (2.31)  evaluated  at  T  is  in  order--cf -magnitude  agree- 

+  ^2 

r.*eno  with,  the  data  reported  by  Atkirs"'  . 
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